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1 INTRODUCTION 

• I— I 

: HU Aqr is a member of the AM Her sub-class of cataclysmic 
5_j ' variables, also called polars because of their highly polarised 

' optical emission. Polars are characterised by a white dwarf 
primary with a strong magnetic field (~ 10 — 240 MG) and 
no accretion disc (see Cropper 1990 and Warner 1995 for 
reviews). The secondary is a main sequence star which fills 
its Roche lobe, resulting in an overfiow of material from 
the inner Lagrangian point to form an accretion stream. 
The stream is initially on a ballistic freefall trajectory, until 
at some radius the magnetic pressure of the white dwarf 
magnetic field is approximately equal to the ram pressure 
of the stream. At this point the stream is threaded onto the 
field lines of the white dwarf and accretes onto a small area 
near one or both of the magnetic poles. 

HU Aqr is an eclipsing polar with a period of 
~ 125 mins, lying just below the period gap (see e.g. Howell, 
Nelson & Rappaport 2001). Harrop-AUin, Hakala & Cropper 
(1999b) used the eclipse of the white dwarf primary by the 
secondary to establish the distribution of brightness along 
the stream given a number of parameters defining both the 
geometry of the system and the physical parameters of the 
model. The stream emission was calculated using model hght 
curves which are optimized by a genetic algorithm, based on 



ABSTRACT 

We apply an eclipse mapping technique to observations of the eclipsing mag- 
netic cataclysmic variable HU Aqr. The observations were made with the S-Cam2 
Superconducting Tunnel Junction detector at the WHT in October 2000, providing 
high signal-to-noise observations with simultaneous spectral and temporal resolution. 
HU Aqr was in a bright (high accretion) state {V = 14.7) and the stream contributes 
as much to the overall system brightness as the accretion region on the white dwarf. 
The stream is modelled assuming accretion is occuring onto only one pole of the white 
dwarf. We find enhanced brightness towards the accretion region from irradiation and 
interpret enhanced brightness in the threading region, where the ballistic stream is 
redirected to follow the magnetic field lines of the white dwarf, as magnetic heating 
from the stream— field interaction, which is consistent with recent theoretical results. 
Changes in the stream eclipse profile over one orbital period indicate that the magnetic 
heating process is unstable. 

Key words: accretion, accretion discs — binaries: eclipsing — novae, cataclysmic 
variables — stars: individual: HU Aqr — stars: magnetic fields 



a method first employed on photometric eclipse profiles of 
HU Aqr by Hakala (1995). The technique was developed and 
improved by Harrop-AUin et al. (1999a) and applied to ob- 
servations of HU Aqr in both a high and low accretion state 
(high- and low-mass transfer; Harrop-AUin et al. 1999b and 
Harrop-AUin, Potter & Cropper 2001). 

This method of eclipse mapping has also been applied 
to emission-lines (Sohl & Wynn 1999; Vrielmann & Schwope 
2001; Kube, Gansicke and Beurmann 2000). Kube et al. 
(2000) used HST FOS spectra of (another polar) UZ For 
to map the accretion stream line emission at Civ A 1550A, 
and Vrielmann & Schwope (2001) used the emission line 
light curves of H/?, H7 and He 11 A4686A of HU Aqr. Both 
these line-emission methods use a three-dimensional stream 
so data from the entire orbit of the system is fitted and the 
model can reproduce the out-of-eclipse features, such as the 
pre-eclipse dips (Watson 1995). They can also reproduce the 
effect of differing brightness distributions between the irra- 
diated and un-irradiated faces of the stream, and this effect 
is not reproduced in the model of Harrop-AUin et al. (1999a 
and 1999b). However, the Harrop-AUin et al. technique has 
the benefit of being sensitive to the total emission from the 
stream (not just the line emission). The technique also has 
the advantage of access to high signal-to-noise ratio data, 
which is crucial in producing robust fits. 
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Figure 1. White liglit curves and sky-backgrounds for cycle 
29982 (upper plot), 29994 (middle plot) and 29995 (lower plot). 
The background is taken from pixels (1,1) and (6,6), while the 
light curves are from all 36 pixels. The light curves are in 1 s 
time bins (0.00013 phase) and the backgrounds in 18 s time bins 
to facilitate direct comparison. The abscissa scales have been ad- 
justed so the eclipses lie at the same locations. A neutral density 
filter was used for cycle 29994, hence the lower count rate and 
increased noise, and the seeing was poor for the first part of cycle 
29982. 



The light curves presented here were obtained using the 
S-Cam2 Superconducting Tunnel Junction (STJ) camera, 
developed by the ESA Astrophysics Division at ESTEC. The 
camera is the second prototype of a new generation of detec- 
tors that record the energy as well as the position and time 
of arrival (to within 5 /xs for this particular detector) of 
the incident photons (Rando et al. 2000). The application of 
STJs to optical photon counting was first proposed by Fer- 
ryman, Foden & Feacock (1993), and has since been applied 
to observations of the Crab Fulsar (Ferryman et al. 1999), 
the magnetic cataclysmic variable UZ For (Ferryman et al. 
2001) and to quasar spectroscopy (de Bruijne et al. 2002). 

The detector itself consists of a liquid helium cooled ar- 
ray of 6 X 6 pixels, each being 25 x 25 ^rr? , corresponding 
to 0.6 X 0.6 arcsec'^ per pixel and a field-of-view of about 
4x4 arcsec^. The pixels are sandwiches of superconduct- 
ing tantalum with a thin insulating layer in the middle and 
the whole device is cooled well below the superconductor's 
critical temperature (about O.lTc). An incident photon then 
perturbs the device equilibrium and as the energy gap be- 
tween the ground state and the excited state is only a few 
meV, a large number of free electrons is created by each pho- 
ton, this number being proportional to the photon energy. 
This is in contrast to normal optical CCD semiconductors 
where the band gap is ~ 1 eV, and photon absorption results 
in typically only one free electron being created. 

The S-Cam2 instrument is particularly suited to obser- 
vations of cataclysmic variables. The high-time resolution 
and simultaneous observations of spectral and intensity vari- 
ations is ideal for eclipsing systems with orbital periods of 
the order of those in cataclysmic variables. The intensity 
variations over the rapid ingress and egress can be probed 
directly, and the rapid variations of the intensity of the ac- 
cretion stream can be used to provide information on the 
possible mechanisms of emission along the stream. 



2 OBSERVATIONS 

S-Cam2 was mounted at the Nasmyth focus of the William 
Herschel Telescope during 2000 October 2/3 and 3/4. Ta- 
ble |l| gives the cycle numbers (relative to the ephemeris of 
Schwope et al. 2001), start time in UTC, time of mid egress 
and the observation length in seconds. 

A total of five eclipses of HU Aqr were recorded. 
The ingress was missed for eclipse 29982 (ephemeris of 
Schwope et al. 2001), and only just recorded for 29993. The 
relatively high count rate of HU Aqr caused some difficulties 
for the data acquisition system, limiting the duration of the 
runs. The first three eclipses exceeded the data acquisition 
system limits, resulting in loss of absolute time reference. 
Eclipse 29994 was therefore taken with a neutral density fil- 
ter with a throughput reduction factor of 10. This is the most 
complete of the eclipses, but suffers from reduced signal-to- 
noise ratio. The data gap in cycle 29982 was caused by the 
problems with the data acquisition system noted above. The 
seeing was in the range 1 to 1.5 arcsec, except at the begin- 
ning of the run for cycle 29982 and throughout cycle 29983, 
when it was 2 arcsec, which led to some spillage of light from 
the 6x6 array. We have therefore not used those data. 

Observations of the spectrophotometric flux standard 
BD+28 4211 were taken to calibrate the data. They were 
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Table 1. Summary of observations of HU Aqr, where cycle number refers to the ephemeris of Schwope et al. (2001). 



Cycle 


Date 


Start time 


Observation 


Egress time 


Remarks 


number 


UTC 


UTC 


length (s) 


BJD (TDB) 




29982 


2000 Oct 2 


21:11 


1420 




Data gap; poor seeing 


29983 


2000 Oct 2 


23:11 


1563 




Poor seeing; data not used 


29993 


2000 Oct 3 


20:13 


831 




Truncated ingress 


29994 


2000 Oct 3 


22:02:42 


2400 


2451821.441021 


NDl filter 


29995 


2000 Oct 4 


00:17:00 


1216 


2451821.527841 


Truncated egress 


BD+28 4211 


2000 Oct 2 


23:46 






Standard; ND2 filter; poor seeing 



made through a neutral density filter with attenuation factor 
100 and also sufi^ered from poor seeing. 



3 DATA REDUCTION 

The data reduction process employs specific pipeline pro- 
cessing of the S-Cam2 data (de Bruijne et al. 2001), based 
on the FTOOLS suite of software (Blackburn 1995) and a full 
description is given by Ferryman et al. (2001) for their ob- 
servations of UZ For. 

Once reduced, the data can be split a posteriori into dif- 
ferent energy bands (wavelength ranges). The intrinsic STJ 
resolution is such that A/AA ~ 9, but here we use only three 
bands to represent 'red', 'yellow' and 'blue' light. The energy 
bands are selected so that roughly equal numbers of events 
occur in each. The wavelength range is limited by the at- 
mosphere at short wavelengths and the use of optics to sup- 
press infrared photons at long wavelengths, and corresponds 
to approximately 340—680 nm. The wavelength ranges for 
the blue, yellow and red bands here are 340—470 nm, 470- 
—550 nm and 550—680 nm. 

As part of the pipeline processing, the data were fiat- 
field corrected using a single map derived from sky observa- 
tions. The data were then corrected for atmospheric extinc- 
tion using the standard La Palma tables of extinction values 
as a function of wavelength and air mass, and finally sky- 
background subtracted. The pipeline processing is applied 
to a fight curve in arbitrary time bins of 1 s. The extraction 
of the background-subtracted light curves used the entire 
6x6 array for the object and the two corner pixels (1,1) 
and (6,6) for the background. The background light level 
is then taken as the mean level during eclipse and this is 
subtracted from the source light curve. 

Figure Q shows the white light curves for eclipses 29982, 
29994 and 29995 before subtraction of the background, to- 
gether with the backgrounds taken from the two corner pix- 
els (1,1) and (6,6). 

Once the data had been calibrated and reduced, we 
folded the data on the orbital period. The linear ephemeris 
of Schwope et al. (2001) was used, which defines inferior con- 
junction of the secondary as (j) — 1.0. The UTC times are 
transformed to TDB (at the solar system barycentre, i.e. 
including light travel times). For cycles 29982 and 29993, 
where there is no absolute time reference, the egress is 
aligned to be at the same phase as in cycles 29994 and 29995. 

Figures ^ and ^ show the white, red, yellow and 
blue light curves, with the colour ratios yellow/red and 



blue/yellow, for cycles 29993, 29994 and 29995. The sec- 
ondary star has been subtracted from the colour ratios (but 
not the light curves) so that the change in the ratios is due 
wholly to the stream, after accretion region ingress. 

We extracted intervals of good seeing from the standard 
star observation to determine the countrate in the yellow 
band, corresponding most closely to the V band. This gave 
a zero point magnitude for 1 count /s in this band of 24.0, 
taking into account the (assumed) factor of 100 from the 
neutral density filter. The corresponding maximum bright- 
ness at (/) = 1.1 in Figure ^ is V = 14.7. This indicates that 
HU Aqr was in a high accretion state at the time of these 
observations (see Schwope et al. 2001). 

4 THE LIGHT CURVES 

The light curves in Figures I and I show a number of features 
that are characteristic of an eclipsing polar system. The most 
prominent feature is the eclipse itself, which starts with the 
limb of the secondary star eclipsing the accretion stream. 
At about 0=0.964 the bright accretion region on the white 
dwarf is eclipsed in a few seconds (the white dwarf is also 
eclipsed around these phases, but is much fainter). After this 
only the secondary star is visible, and the sequence is then 
approximately reversed on egress. 

4.1 Pre-eclipse dip 

Prior to the eclipse of the white dwarf, we observe a pre- 
eclipse dip. Figure ^ shows this dip centred at <jf> « 0.872. 
The dip is caused by the eclipse of the accretion region on 
the white dwarf by the magnetically confined section of the 
accretion stream. Consequently the phase of the centre of the 
dip is directly related to the azimuth of the coupling region, 
where the stream is threaded onto the field lines of the white 
dwarf. The centre of the dip in eclipse 29994 corresponds to 
an azimuthal angle of 46'' which is within the range of values 
found by Harrop-AUin et al. (1999b) of 42° to 48° for their 
light curves. 

Schwope et al. (2001) investigated the correlation be- 
tween the phase of the dip and the soft X-ray luminosity, and 
concluded that the dip ingress occurs earlier in phase when 
the system is brighter. The ingress of the dip at (/> « 0.856 
indicates that the soft X-ray luminosity is at the higher end 
of the range, as expected for the bright state of the system. 

The dip is visible only in cycle 29994, the other observa- 
tions being too short to cover the relevant phases. We can- 
not therefore relate its movement between successive eclipses 
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Figure 2. The sky-subtracted phase folded light curves for eclipses 29993 (left) and 29995 (right), representing from the top white, red, 
yellow and blue light, with a resolution of 1 s, or 0.00013 phase. The colour ratios are binned in 5 s bins and the contribution from the 
secondary during total eclipse has been removed. Points with errors over 0.5 in the colour ratio plots have been omitted for clarity - 
these occur during complete eclipse when there is effectively zero count rate after subtraction of the secondary contribution. To facilitate 
comparison, dotted lines have been placed at the arbitrary levels of 0.8 and 1.2 in the blue/yellow and yellow/red ratios respectively. 
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Figure 3. As for Figure g but for eclipse 29994. The countrate is lower by a factor of 10 because of the neutral density filter. 
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to t he c hange in stream ecUpse profile as discussed in Sec- 



tion 4.3 



4.2 Accretion region eclipse 

The soft X-ray data in Schwope et al. (2001) suggest that 
there is only one accreting pole. It is possible that there is 
a second accreting pole, as discussed in Harrop-AUin et al. 
(1999b). However, there is no evidence for such a pole in 
our eclipse profiles, so for the purpose of our modelling we 
assume accretion onto one pole. 

Around (j> ~ 0.964 the accretion region has been com- 
pletely eclipsed and the accretion stream is thereafter the 
dominant source of the observed brightness with a small 
contribution in the red band from the secondary. We show 
in Figure ^the accretion spot profile on ingress and egress. 
The profile is constructed by subtracting successive intensi- 
ties over 1 s time intervals (0.00013 phase) and so shows the 
rate of change of the light curves intensity. The profile on 
egress is the mean of cycles 29982, 29993 and 29995 and is 
therefore more reliable than the ingress profile, which uses 
only cycle 29995. The profile may be asymmetrical, with the 
leading edge (corresponding to later phases) of the cyclotron 
emitting region being brighter than the trailing edge. The 
accretion spot ingresses last for ~ 7 s and the egresses are 
~ 8 s. This is unchanged from the durations measured by 
Case (1996) in the data used by Harrop-AUin et al. (1999b) 
and should be compared to the duration of 1.3 s for the soft 
X-ray eclipses (Schwope et al. 2001), so that the extent of 
the cyclotron emitting region is ~ 15° in the optical com- 
pared to ~ 3° in soft X-rays. From soft X-ray absorption 
modelling, Schwope et al. (2001) suggest that the threading 
region is extended in azimuth, but the soft X-ray emission 
originates over only a small range of azimuths: the larger re- 
gion seen in the optical is therefore likely to be the result of 
this extension in azimuthal threading. The cyclotron emit- 
ting region is significantly more extended than that mea- 
sured for UZ For, where Ferryman et al. (2001) found both 
regions extended over 3° ± 1°. 

In those cycles where we have absolute time infor- 
mation (29994 and 29995) the egress takes place at = 
1.0413 in both cycles, 10 s later than predicted by the 
linear ephemeris. This indicates that the observed-minus- 
calculated (O-C) residuals in Schwope et al. (2001) Figure 
4 are continuing to increase. In case there are systematic 
differences between the optical and soft X-ray egress times, 
we rephased the Harrop-AUin (1999) data in their Figure 1 
on the new ephemeris, including the UTC-TDB time cor- 
rection of 59 s appropriate to the Schwope et al. (1997) 
ephemeris epoch. The data were taken on cycles ~ 1322 
in Figure 4 of Schwope et al. (2001), and the O-C of -1.6 s 
is in agreement to within ~ 1 s with those determined from 
the soft X-ray data. The continuing increase in O-C residu- 
als is an indication that the period in the linear ephemeris of 
Schwope et al. (2001) is slightly too short, and is in the oppo- 
site sense to that expected from their quadratic ephemeris 
which predicts a decrease in the O-C residuals. A period 
of 0.0868204212 days would eliminate the 10 s residuals we 
observe. 

The duration of the eclipse from the centre of the ac- 
cretion region ingress to the centre of its egress is = 
0.0778 ± 0.0002 in cycles 29993, 29994 and 29995. This 






Figure 4. (Top) The accretion spot profile on ingress (from cycle 
29995) and egress (the average of cycles 29982, 29993 and 29995). 
The profiles are constructed by subtracting successive intensities 
over 1 s time intervals. The ordinate of the upper plot is inverted 
to facilitate comparison. (Bottom) The view of the eclipsing limb 
at accretion region mid-ingress and mid-egress. The accretion re- 
gion subtending 15° is drawn to scale. The radius of the white 
dwarf is assumed to bo 0.017a, where a is the binary separation. 



is the same as that found by Harrop-AUin et al. (1999b) 
of = 0.0779 ± 0.0002, and Schwope et al. (1997) of 
IS.(j) = 0.0779. The lack of variation in the duration of the 
eclipse is remarkable and implies that the accretion spot is 
at the same latitude for these observations. On the other 
hand, this is shorter by 6 s than the 0.0782 - 0.0792 mod- 
elled by Schwope et al. (2001) from the soft X-ray eclipses, 
where the ingress was uncertain due to the low countrate 
and absorption. 

For a particular mass ratio, the width of the eclipse de- 
termines the inclination of the system. Using q — M2/M1 = 
0.25 (Schwope et al. 2001), the inclination is 85.0° (assum- 
ing Ml = 0.9Mq and /3 = 40°, C = 50°, where Mi is the 
mass of the white dwarf and /3, ^ are the magnetic longitude 
and colatitude of the accretion region. These correct for the 
location of the accretion spot on the surface of the primary; 
see below). 



4.3 Stream eclipse variations 

Complete ingress of the stream occurs at = 
0.982, 0.979, 0.977 and 0.982 for eclipses 29982, 29993, 29994 
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Figure 5. The ingress and egress of eclipses 29982 (dotted), 29993 
(dash-dot), 29994 (soUd) and 29995 (dashed) smoothed over a 
running average of 20 s. The ordinates are normahsed to 1 at the 
end of spot egress, and to at total eclipse. 

and 29995 respectively. During the period of complete 
eclipse the secondary is the only contributor to the hght 
curve and provides a constant contribution which is great- 
est in the red band, as expected, with no contribution in the 
blue band. 

Figure ^ shows more clearly that the observations reach 
total stream eclipse at different phases even though the spot 
ingresses/egresses both occur at the same time. This rapid 
change in the accretion flow— magnetic field interaction be- 
tween successive eclipses was first observed by Glenn et al. 
(1994) who found a difference in the time it took for the 
stream to be completely eclipsed of more than a minute be- 
tween their two successive eclipses. 

The shape of the accretion stream eclipse profile in the 
observed hght curves gives a qualitative prediction of where 
brightness enhancements can be expected in the accretion 
stream. When the stream is brighter at later phases in the 
stream eclipse, as in cycles 29982 and 29995 (Figure^), there 
must be more emission in the threading region, as this is the 
only part of the stream still visible. For the same reason, 
if the final stage of the stream eclipse takes place at later 
phases, the threading region must have moved: thus in cy- 
cles 29982 and 29995, either the accretion stream is brought 
further out from the line of centres during the ballistic part 
of the trajectory, or the magnetically channeled stream rises 
further out of the orbital plane (we address this in Section ^ 
below). This could be caused either by a small change in 
the location of the accretion spot on the primary surface, 
moving the accretion spot further from the line of centres, 
or perhaps a change in the velocity of the stream as it leaves 
the Li point. Likewise, if the accretion stream is fainter near 
the white dwarf, the system brightness will be fainter imme- 
diately after accretion region ingress. This may be the case 
for cycle 29995 (Figure |). 

4.4 Colour variations 

For cycle 29995 there is a rise in the yellow/red ratio dur- 
ing the stream ingress (Figure |^, so that successively more 
red light is blocked from the stream as the eclipse pro- 



gresses. This implies that the last part of the stream to be 
eclipsed, near the threading region, emits relatively less at 
longer wavelengths than the rest of the stream. This rise 
is not evident in cycle 29993, but maybe present at a re- 
duced level in cycle 29994 (Figure ^ . This trend is consis- 
tent with a threading region which is increasingly hotter in 
the three successive eclipses. At the same time as the large 
yellow/red ratio change in cycle 29995 there is no change in 
the blue/yellow ratio. 

After eclipse when the stream close to the primary is un- 
covered first, the different yellow/red ratios in cycles 29993 
and 29995 are consistent with a hotter stream near the 
white dwarf in eclipse 29993 than in 29995. An increase in 
the blue/yellow ratio after eclipse is seen in all three cycles 
(29993, 29994 and 29995) and also implies that the magnet- 
ically confined stream material is becoming hotter towards 
the threading region. 

The pre-eclipse dip around « 0.87 in cycle 29994 
(Figure H) is also slightly bluer than at earlier or later 
phases. This was also the case in the hght curves of Harrop- 
AUin et al. (1999b). 



5 INDIRECT IMAGING 
5.1 Methodology 

We have inferred from a qualitative analysis of the light 
curves that the trajectory and/or intensity of the stream 
varies over successive orbital periods, and we can identify 
brightness enhancements along the stream length. However, 
it is more difficult to infer the location of the brightness 
enhancements without a model for the stream trajectory. 
In particular, at stream ingress there is no information on 
which parts of the stream are eclipsed at which phases. We 
can use an eclipse mapping method to create a model accre- 
tion stream by placing model stream points along a prede- 
termined trajectory, consisting of a ballistic part from the 
Li point, coupling to a magnetically confined section where 
the material is threaded by the field lines of the white dwarf. 
The distance from the white dwarf where this transition oc- 
curs is the threading radius R/j,, and is a user input to the 
model. The model system is rotated while the Roche lobe 
filling secondary eclipses the stream and white dwarf, with 
the brightness of the visible model stream points at each 
phase summed to form a model light curve. The resultant 
light curve is optimized with the genetic algorithm (GA) 
evolving the best fit light curve. 

The 'goodness of fit' of a model light curve is measured 
with a 'fitness function' (see Ifarrop-AUin et al. 1999a for a 
full description) . The GA adjusts the brightness points along 
the stream in an attempt to minimise this function, which 
consists of a term and a maximum entropy term, which 
ensures the problem is not under-constrained. After the fi- 
nal stages of the GA, a more conventional line-minimisation 
routine (Powell's method) is used to reach the final mini- 
mum of the solution. 

The accretion stream itself is taken to be physically thin 
in that it does not eclipse the primary, and so features such 
as the pre-eclipse dip (e.g. Watson 1995) are not reproduced 
by the model. However, the brightness contribution of each 
stream point is taken as the sine of the angle between the 
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line of sight and the tangent to the stream at that point, 
essentiaUy a projection effect, thus optically thick. We also 
exclude the white dwarf from the model, although the egress 
of the white dwarf can be seen in the light curves immedi- 
ately before the egress of the accretion region. 

The exact method we apply for the modelling is dif- 
ferent to that of Harrop-AUin et al. (1999b, 2001) in that 
we are not applying the technique to the complete eclipse 
light curves, nor can we interpret the results in the same 
manner. There are a number of reasons for this. Firstly, and 
importantly, our observations are all truncated, either in the 
ingress or the egress of the accretion stream. For observa- 
tions where the egress is truncated (cycle 29995) , the model 
has difficulty breaking the ambiguity which exists when as- 
signing brightness to points eclipsed in the same phase in- 
terval, i.e. the brightness can be either in the magnetically 
confined region or the ballistic. For observations where the 
pre-eclipse light curve is missing, we lack information on the 
accretion region, and parts of the stream towards the sec- 
ondary (cycle 29993). A further difficulty lies in determining 
the nature of the variable emission from the accretion region. 
This is particularly important at egress, where the bright- 
ness of the spot can greatly affect the brightness distribution 
towards the white dwarf. After the egress of the accretion 
region the light curve consists of a variable component, plus 
the stream, which changes rapidly over the course of the ob- 
servation. We have therefore chosen to use a restricted mod- 
elling technique, which relies on the ingress of the accretion 
stream alone. This means that we have truncated the light 
curves, removing all phases up to and including the ingress 
of the accretion region, and those after = 1.0. This is the 
only part of the light curve unaffected by the variable emis- 
sion from the accretion region. We still have the problem of 
the ambiguity of the stream points, but we can make some 
progress in interpreting our results by referring back to the 
light curves and the colour ratios. 

5.2 Fixed parameters and 7?^ 

A number of parameters is required to produce a model 
light curve and reproduce the stream brightness distribu- 
tion. These fall into two types: physical parameters such as 
the masses of the two component stars, and geometric pa- 
rameters such as the location of the accretion region on the 
white dwarf (see Harrop-AUin et al. 1999a for details). The 
model is particularly sensitive to the exact value chosen for 
the parameter i?^, but we can constrain the value using the 
model stream geometry. If the value is too large then the 
model needs to assign a large amount of brightness to a 
few points. Conversely if is too small then no emission 
is assigned by the model to the points in the threading re- 
gion. From brightness maps for a range of values of we 
can determine the end point of the ballistic trajectory and so 
provide a constraint on the value of . We have found from 
our model fits of the light curves used here that the tech- 
nique depends heavily on the data being of a sufficiently high 
signal-to-noise ratio, which is important because of the sen- 
sitivity to the value of R^. Therefore, in order to determine 
the best fit value to use we further restrict the application of 
our model technique to cycles 29993 and 29995 (Figure ^). 

Figure ^ shows model fits to the 'blue' light curve of cy- 
cle 29995 for different values of _R^. The stream maps high- 



Table 2. The values of the coupling radius (-R^i) and location of 
the accretion spot (/3 and f , magnetic longitude and colatitude 
respectively) for cycles 29993 and 29995. Values of f) and have 
an estimated range of ±10°(1(t). 



Cycle number 






c 


29993 


0.26a 


55° 


30° 


29995 


0.22a 


60° 


40° 



light the dependence of the fits on a correct value for 
For a value of = 0.18a the brightness of the threading 
region is low. This lack of emission, or 'hole', is caused by i?^ 
being too small - the data are incompatible with emission 
at the end part of the resulting long ballistic stream. On the 
other hand, if R^ is too large, so that the model ballistic 
stream is shorter than it is in reality, a pile up of excess 
brightness at the threading region is seen, as in the stream 
map for i?^ = 0.26. The inappropriateness of the latter can 
be deduced from the poor fit to the region at = 0.98 (cor- 
responding to the threading region) for values of R^,, which 
are too large (insets to Figure h) . 

For all our modelling we used fixed values of M\ = 
O.9M0, q = 0.25 and i = 85.0° (Harrop-AUin et al. 1999b, 
Schwope et al. 2001). Table ^ gives the values of R^ and the 
field orientation parameters for each cycle. We estimate the 
range on the values of /3 and ^ as ±10°(lcr). 

The values of the three parameters in Table ^ are dif- 
ferent from those found by Harrop-AUin et al. (1999b). We 
find Rf, = 0.22 - 0.26a (1.3 x 10^° cm to 1.6 x 10^° cm), 
whereas Harrop-AUin et al. (1999b) found = 0.17-0. 20a. 
We can attribute this difference to an improved stream tra- 
jectory compared to that of Harrop-AUin et al. (1999a). 
This arises from an incomplete treatment of the appropriate 
forces in the frame of reference chosen by Harrop-AUin et al. 
(1999a). The difference is small, but for the angles used here 
the resulting R^ is significantly different, with the Harrop- 
AUin et al. trajectory underestimating the values for i?^. 

The decrease in the value of i?^ from cycle 29993 to 
29995 implies that the stream penetrates further into the 
magnetosphere before threading onto the magnetic field lines 
(Figure ^) . The effect of the change in the threading radius 
can be seen in Figure ^. The values of /3 ~ 60° , C, ~ 40° found 
here are also larger than those used by Harrop-AUin et al. 
(1999b), but consistent within the uncertainties. Any change 
in these parameters would imply that the location of the 
accretion spot has moved further from the line of centres 
and the altered geometry of the magnetic field carries the 
accretion stream further round the white dwarf from the line 
of centres (see Section 5). They are in agreement with those 
found by Schwope et al. (2001) (their figure 6). 



5.3 Model results 

With the values of i?^ determined above, we can use the 
model to find the best fit brightness map, from which we 
can then determine the brightness per unit stream length 
along the stream trajectory. This will show the brightness 
independent of the actual length of the stream eclipsed in 
any particular phase interval. This is important as there 
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Figure 6. Fits to the 'blue' light curve of cycle 29995 with different values for the threading radius (0.18a, 0.22a and 0.26a). The 
stream brightness maps on the right are gray scale plots with brighter stream emission being represented as darker. The gradual closing 
of a 'hole' at illustrates the sensitivity of the model to its location. The insets in the left hand panels show a close-up view of the fit 
at the threading region. For this cycle the value of i?^ = 0.22a represents the best stream trajectory. 



are difFerent lengths of ballistic and magnetically confined 
stream eclipsed in a given phase interval. 

Figure |^ presents the results of the modelling in two 
columns representing the two cycles presented here and the 
three rows the three energy bands: red, yellow and blue. We 
have divided the stream into three sections, as illustrated in 
Figure ^, which enclose the section nearest the white dwarf, 
the threading region and the stream in between the two. 
The ordinate then represents the total brightness of points 
eclipsed in these intervals per unit stream length, that is 
points in the ballistic and magnetically confined sections of 
the stream eclipsed in that phase range. The results have 
been calibrated into energy units using the standard star 
observation in order to facilitate comparison between the 
different bands. However the calibration is only approximate 
due to a lack of a standard star on the second night of ob- 
servations. 

We can see from Figure ^ that there is enhanced bright- 
ness in the section containing the threading region in all 
bands in cycle 29993, but only in the yellow and blue in cy- 
cle 29995. The section containing those parts of the stream 
nearest the white dwarf are also bright in cycle 29993. How- 
ever, in cycle 29995 this is only seen in the blue. The bright- 



ness of the difFeren t se ctions is consistent with the colour 
ratios from Section 4.4. In Figure |^ the threading region in 
cycle 29995 is hotter than that in cycle 29993 as the ratio 
of the yellow to red is greater for cycle 29995. This is again 
consistent with the colour ratios (Section 4.4). However, the 
blue to yellow ratio is similar in the two eclipse ingresses 
suggesting that the increase has occurred in both blue and 
yellow bands. 



6 DISCUSSION 

6.1 Comparison with earlier results 

The model fits of Harrop-AUin et al. (1999b) and Harrop- 
AUin (1999) for one pole accretion, show a general en- 
hancement in the threading region and towards the white 
dwarf. Some eclipses (for example cycle 3723 and 3724 in 
Harrop-AUin 1999) show significantly greater brightness in 
the threading region of the white dwarf than others (such as 
the immediately preceding cycle 3722). 

Enhanced brightness regions are also found by Vriel- 
mann & Schwope (2001) and Kube et al. (2000) who apply 
similar modelling techniques. Vrielmann & Schwope (2001) 
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Figure 7. The above plot shows for cycles 29993 and 29995, in 
each band, the total brightness per unit stream length in three 
sections of the stream. These intervals include a section of the 
ballistic trajectory plus those parts of the stream nearest the 
white dwarf (section 1), the threading region (section 3) and the 
magnetically confined stream between the two (section 2), c.f. 
Figure ^. A brightness enhancement at the threading region is 
evident in both cycles, but an enhancement towards the white 
dwarf is not seen in the red and yellow bands of cycle 29995. 




Figure 8. The upper plot shows the model stream points for 
cycles 29993 and 29995 looking parallel to the orbital plane, while 
the lower plot shows them looking down onto the orbital plane. 
The difference in the stream trajectories is evident. The shaded 
regions are three equal phase intervals encompassed by the ingress 
of the white dwarf and the final ingress of the threading region 
for cycle 29993. These correspond to those parts of the stream 
used to form sections 1, 2 and 3 in Figure ^ (For cycle 29995 the 
definitions are equivalent.) 



apply their Accretion Stream Mapping technique to H/3, H7 
and He 11 A 4686A observations of HU Aqr. They find a 
brightness enhancement at the threading region, but not to- 
wards the white dwarf, similar to cycle 29995, so implying 
that this may be absent in their emission lines. Kube et al. 
(2000) used the Civ A 1550A line emission from UZ For, with 
a 3-dimensional stream model. They found three regions of 
enhanced brightness: one on the ballistic accretion flow, and 
two on the magnetically confined section. They suggest that 
the enhancements on the magnetically confined section are 



caused by irradiation of denser sections of the stream with 
a large area near the accretion region, and a smaller region 
near to the threading region. This may be the case in cycle 
29993, where we find a bright stream near the white dwarf, 
and again near the threading region. Although Kube et al. 
find no enhancement actually at the threading region, this 
may be a result of an increase in the density of the stream 
as it approaches this point, resulting in an increase in the 
continuum optical depth, and hence a decrease of the Civ A 
1550 A equivalent width, and so is not necessarily indicative 
of a faint stream at this region. 



6.2 Heating in the magnetically confined stream 

Heating of the magnetic section of the stream has been mod- 
elled by Ferrario & Wehrse (1999) for a stream which is as- 
sumed to thread onto the field lines at a coupling radius Vc 
(our Rfj) from the white dwarf, over a radial distance Arc in 
the orbital plane. This provides an opportunity for a com- 
parison between our observationally derived results, and the 
main points of their theoretical model results. 

Ferrario & Wehrse (1999) consider two heating mech- 
anisms: irradiation by the X-ray component from the ac- 
cretion shock (the soft X-rays being up to ~ 4000 times 
more efficient at heating the stream than the hard X-rays) 
and the effects of magnetic reconnection in the stream-field 
interaction over the threading region Arc. The comparison 
we make with these theoretical models is qualitative as the 
geometrical structure of the magnetically confined accretion 
fiow in Ferrario & Wehrse is funnel-shaped, compared to our 
linear trajectory. However, our results showing brightness 
enhancements towards the threading region are consistent 
with their models which incorporated magnetic heating in 
the threading region. This implies that some magnetic heat- 
ing mechanism is needed. 



6.3 Temporal variations in the stream profile 

Based on the colour ratios we deduced that the threading 
region was brig hte r and hotter in cycle 29995 than in cycle 

29994 (Section 4.4). There are therefore significant changes 
in the threading region on the timescale of the orbital period 
(125 mins). 

Although we have only modelled two cycles here, we can 
see that there is a difference in the brightness confirming the 
variability seen in the colour ratios and in the stream eclipse 
profiles from the light curves (see Figure ^ . The large bright- 
ness enhancement near the white dwarf in cycle 29993 could 
be irradiated material which cools over the timescale of the 
next orbital period, leaving a fainter region as seen in cycle 

29995 where the observed enhancement is less pronounced. 
However, without the model of the intermediate cycle, which 
we have excluded for reasons discussed previously, we cannot 
be certain and it will require more high signal-to-noise ratio 
observations of consecutive eclipses to investigate whether 
this in fact occurs. What is clear from the model results, 
and indeed directly from the raw light curves and colours, is 
that the emission from the whole stream is highly dynamic 
and unstable. 

Dramatic changes in the stream eclipse profiles were 
also observed in the low state (Harrop-Allin et al. 2001, 



© 0000 RAS, MNRAS 000, 000-000 



ST J Observations of HU Aqr 11 



Glenn et al. 1994). These cycle-to-cycle changes in stream 
brightness and trajectory require a revision in our view of 
the stream, and in the manner in which the magnetic heat- 
ing takes place in the threading region. Future treatments 
may need to consider large scale magnetic instabilities, and 
quasi-cyclic behavior. 
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7 SUMMARY AND CONCLUSIONS 

We have carried out high signal-to-noise ratio observations 
of HU Aqr using S-Cam2 on the WHT on two nights. The 
system was in a high accretion state, and from the single 
sharp change in the eclipse profile and archive soft X-ray 
light curves in this state, we infer that matter was accret- 
ing at only one pole on the white dwarf. At the onset of 
the eclipse, the accretion stream is the source of more than 
half of the optical emission from the system. The system 
brightness was similar from orbit to orbit, and the eclipse 
duration remained constant, but the shape of the accretion 
stream eclipse changed significantly, ending ai (p — 0.979, 
0.977 and 0.982. 

We find eclipse durations which are unchanged from 
past optical studies (Harrop-AUin et al. 1999b, Schwope 
et al. 1997), but shorter than those deduced in the soft X- 
rays by Schwope et al. (2001). However, the location of the 
accretion region is similar to that found by Schwope et al. 
(2001), so that there is no evidence that the optical emission 
is from higher latitudes on the white dwarf than the soft 
X-rays. This indicates that the duration of the soft X-ray 
eclipses may be affected by absorption, as they suggested. 

The duration of the egress of the accretion region in 
the optical is 8 s, compared to 1.3 s in soft X-rays (Schwope 
et al. 2001). This indicates clearly that the region emitting 
cyclotron radiation is extended by a factor of ~ 5 by com- 
parison with the soft X-ray emitting region, which Schwope 
et al. (2001) calculated as subtending an angle of 3°. 

We have found significant changes in the colour of the 
accretion stream from one eclipse to the next. This indicates 
that the threading region is hottest in the last of th e ec lipses 
by comparison with the previous two (see Section 4.4). We 
have modelled the stream using the technique of Hakala 
(1995) and Harrop-AUin et al. (1999a). This finds that most 
of the emission originates from two places, the region close 
to the white dwarf and in the threading region. By com- 
parison with the models of Ferrario & Wehrse (1999) this 
indicates that magnetic heating is required in the thread- 
ing region. The modelling clearly identifies an increase in 
brightness in the threading region for both cycles and an 
enhancement towards the white dwarf for cycle 29993. From 
this change in brightness of the heated regions in the model 
streams, and the varying stream eclipse profiles, we suggest 
that the magnetic heating in the threading region may be 
unstable. The implications of the highly variable stream tra- 
jectory and brightness profile should be recognised in future 
investigations of the stream properties. A systematic study 
of series of consecutive eclipses, with the highest possible 
signal-to-noise ratio, is required to investigate the charac- 
teristics of the magnetic heating and stream instabilities. 
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